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We analyze the spectrum of the 2D S=l/2 frustrated Heisenberg model near the transition from the 
spontaneously dimerized spin-liquid phase into the Neel ordered phase. Two excitation branches: 
the triplet magnon, and the collective singlet mode, both become gapless at the transition point. 
However we find that the length scales associated with these modes are well separated at the quantum 
transition. While in the quantum disordered phase the singlet excitation has finite spectral weight 
and reflects the existence of spontaneous dimer order, near the transition point the size of the singlet 
bound state grows exponentially with the correlation length, and hence the quasiparticle residue is 
exponentially small. Therefore the critical dynamics remains in the 0(3) universality class in spite 
of the four gapless modes. 
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Quantum phase transitions between magnetically or- 
dered and disordered phases can take place at T = 
by varying the exchange interactions which can drive the 
spin-spin correlations from long-range behavior, charac- 
terized by an infinite correlation length, towards a short- 
range regime, typical for disordered phases. An example 
of a quantum model which exhibits such transitions is 
the two-dimensional (2D), 5=1/2 Heisenberg antifer- 
romagnet (HAFM) on a square lattice. While for uni- 
form nearest-neighbor interactions the HAFM has long- 
range Neel order in the ground state with sub-lattice 
magnetization M sw 0.3 Q|, inclusion of additional in- 
teractions, such as dimerization and/or frustration, leads 
to increased quantum fluctuations and ultimately vanish- 
ing of M at a critical coupling. Examples of transitions 
caused by local alternation of the exchange couplings are 
the dimerized HAFM the two-layer HAFM f||§ and 
the C0V4O9 lattice (1/5-th depleted square lattice) Q. 
In these cases the local dimer or plaquette correlations 
eventually win over the long-range Neel order, leading 
to a non-magnetic ground state. Another route towards 
a magnetically disordered ground state is introduction of 
frustrating second- neighbor interactions ( J2), in addition 
to the nearest-neighbor ones (J\) (see Fig.l.). The Neel 
order disappears at (J2/Ji)c ~ 0.4 in this case |7|-[To|l. 

An important issue concerning the quantum transi- 
tions mentioned above is their universality class. It is 
generally accepted that the effective low-energy theory 
for the 2D Heisenberg systems with a collinear (Neel) 
order parameter is the 0(3) non-linear sigma model 
(NLSM) in 2+1 dimensions Q. This field theory con- 
tains a single effective coupling constant g and, at T = 0, 
describes the ordered Neel phase for g < g c . For g > g c 
the NLSM is in a quantum disordered phase with a finite 
correlation length. However the determination of g c and 
the nature of the disordered phase are beyond the field 
theory formulation and depend on the specific details of 
the model. In addition, Berry phases associated with in- 



stanton tunneling between topologically different config- 
uration are present in the NLSM . In one dimension 
the Berry phase effects are known to be important, essen- 
tially leading to the difference between the excitations in 
the integer and half odd-integer spin chains [[[3) . In 2D 
Berry phases are also present but their role is less clear. 
If one neglects these purely quantum effects, the univer- 
sality class of the quantum transitions in the 2D HAFM 
should be the same as that of the classical 0(3) vector 
model in 3D [14] . Quantum Monte Carlo simulations 



performed on the two-layer HAFM |15| and the CaV^Og 
lattice HAFM |l(| confirm with high accuracy that the 
quantum transitions in the above two models are in the 
0(3) universality class. There also has been a report |l7j] 
that the 2D dimerized HAFM exhibits a deviation from 
the 0(3) behavior, which is presumably due to the small 
size lattices studied in the above work. Generally the 
most accurate Monte Carlo results seem to indicate that 
the quantum Berry phase effects are not important, at 
least in the models where the quantum transitions are 
driven by explicit (exchange driven) dimerization. 

On the other hand the J\ — J2 model, which exhibits a 
quantum transition due to frustration, has a much bet- 
ter chance for deviation from the 0(3) universality class. 
The reason is that the Berry phases were shown to be 
relevant and intimately related to the presence of spon- 
taneous dimer (spin-Pcierls) order in the quantum disor- 
dered phase of this model fll|[iij . Within the formalism 
of the large N expansion for the Sp(N) theory (N = 1 be- 
ing the physical limit), Read and Sachdev jll| found two 
divergent length scales at the transition from the quan- 
tum disordered into the Neel ordered phase. The first one 
is the usual correlation length £ which governs the expo- 
nential decay of the spin-spin correlations in the disor- 
dered phase and is inversely proportional to the (triplet) 
magnon gap, £ cx 1/A. The instanton effects however 
lead to the appearance of spontaneous spin-Peierls corre- 
lations and a second gapped singlet mode with a charac- 
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teristic scale £sp (inverse singlet mass). The two scales 
are related via: 

Zsp~t CN , N»l, (1) 

where C = d + 0(1 /N), G\ ~ 1. Since N is large one 
expects £sp 3> £. The presence of two divergent length 
scales at the transition point would naively suggest a 
change in the universality class. However, as argued in 
Ref. ^0|, the fact that £sp is a (large) power of £, caus- 
ing the two scales to be well separated near the critical 
point, is a characteristic feature of a dangerously irrele- 
vant coupling. This means that even though the Berry 
phases are relevant in the disordered phase, ultimately, 
near the critical point, their effect disappears. In partic- 
ular the dimer order parameter D is expected to behave 
as D ~ £gp ~ A CN and thus vanishes very fast as the 
critical point (A — » 0) is approached |^9| . In this sce- 
nario the quantum critical fluctuations of the Neel order 
parameter are decoupled from the singlet mode and con- 
sequently the transition is still of 0(3) type. Notice that 
the above analysis is certainly valid provided the 1/iV ex- 
pansion behaves well, since only then the N — ► oo results 
are relevant to the physical situation N = 1. However 
corrections beyond the N — oo limit have not been sys- 
tematically calculated in the literature, due to the com- 
plex nature of the problem. 

The purpose of the present work is to analyze the struc- 
ture of the excitation spectrum and the scales that ap- 
pear near the quantum critical point in order to test the 
Sp(N) field theory predictions. We work directly with 
the physical spin problem (TV = 1) and the approxima- 
tion scheme that we use is based on a perturbative expan- 
sion around the spontaneously dimerized ground state in 
the quantum disordered phase. First, let us mention that 
the numerical implementation of this expansion via the 
dimer series expansion l2l]-|2^| , as well as the mean-field 
and diagrammatic treatments P2| , confirm the sta- 
bility of the spontaneously dimerized phase for interme- 
diate values of frustration. This means that the large N 
limit captures the essential physics of the problem, even 
though it can not be trusted numerically in regards to the 
exact location of the phase boundaries. However both the 
series expansions and the diagrammatic method are not 
accurate enough to calculate reliably the critical expo- 
nents near the the transition into the Neel phase, since 
the exponents are not expected to vary considerably. For 
example the exponent v governing the vanishing of the 
triplet gap: A ~ (g - g c ) u , is v » 0.71 for the 0(3) 
and v w 0.75 for the 0(4) universality class (w|. Such 
small difference can not be confidently resolved with the 
above methods pa] . This is why we will follow a different 
route, namely we will analyze the possibility of having ad- 
ditional soft modes at the transition (in addition to the 
triplet mode). 

The Hamiltonian of the frustrated Heisenberg antifer- 
romagnet is: 



H = J 1 ^2 Si-Sj + J2 S - S J' ( 2 ) 

NN NNN 

where J\ > is the nearest-neighbor, and J2 > is the 
frustrating next-nearest-neighbor exchange coupling on a 
square lattice (defined as shown in Fig.l). All the spins 
Si = 1/2. In order to find the excitation spectrum of H 
we follow closely the treatment of Ref. |2^] which is briefly 
outlined below. The starting point is grouping the spins 
into dimers (singlets) in the pattern, shown in Fig.l. This 
configuration (which is degenerate with three others, ob- 
tained by translation by one lattice site, rotation by 7r/2, 
and rotation plus translation) was found to be stable in 
the parameter window (fo/Jilci < J2/J1 < {Jij J\)ci 
p2| . Here (Jij J\) c \ ~ 0.38 is the transition point into 
the Neel phase, the neighborhood of which is the region 
we want to analyze. The Hamiltonian can be expressed in 
terms of bosonic operators t\ a ,a = x, y, z creating three 
degenerate triplet excitations from the singlets formed by 
each pair of spins, as shown in Fig.l. The site index i 
now numbers the sites on the dimerized lattice. The ef- 
fective Hamiltonian describing the interactions between 
the triplets is f2§: 

H = H 2 + H 3 + H 4 , (3) 
H2 = J2 \ A Aj^ + *y (4/-ka + h - c -) ) . ( 4 ) 

#3 = Y R ( k l> k 2)e Q! ^7*k 1 a*k 2 ^*k37 + h - c -' ( 5 ) 
1+2=3 

H 4 = ^ [T( k i - ka)(SaS$p-y - SapS-fs) + (6) 
1+2=3+4 

The following definitions are used in Eqs.(||||): = 
Ji - (Ji/2)&. + ( Ji - J 2 )Cfe 8 - MkAk y ,B k = A 1t -J 1 , 
and the matrix elements in the quartic and cubic in- 
teraction terms: 4T(k) = Ji£ kx + 2(Ji + J2)6c„ + 
2J 2 ^kM y , 4R(p,q) = -Ji7 Px - Wflwipy ~ {p ~> <?}, 
where we have defined £fc = cos(fc), 7^ = sin(fc). The 
T and R terms describe the inter-site interactions arising 
from the exchange between the dimers. An additional 
on-site (U) term is also introduced and one must set 
U — > 00. This term reflects the hard-core nature of the 
bosons which follows from the kinematic constraint on 
the Hilbert space t\ a t\p = 0. The constraint is necessary 
in order to ensure that the bosonic Hamiltonian in terms 
of the triplet operators corresponds uniquely to the origi- 
nal spin Hamiltonian (|^) and no unphysical states appear 
in the final result. The sums over k extend over the Bril- 
louin zone of the dimerized lattice, i.e. — 7r ^ k X f ky 7r. 
In this notation the Neel ordering wave-vector ((jr, 7r) of 
the original lattice) is Qaf — (0,7r). 
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The spectrum of Eq.(||) was studied in Ref. |22j by 
summing selected infinite series in the perturbative ex- 
pansion. The dilute Bose gas approximation was used 
and the diagrams classified in powers of the density of 
magnons. The diagrammatic treatment was also com- 
pared with numerical results obtained by high-order 
dimer series expansions and the agreement was found to 
be very good. We will therefore present only diagram- 
matic results from now on. In the quantum disordered 
phase {J2I J\)c\ < J2/J1 < {J2I J\)c2 the triplet exci- 
tation spectrum w(k) has a non-zero gap A = lo(Qaf) 
which reflects the fact that the dimer configuration is sta- 
ble. As the critical point {J2I Ji)ci = 0.38 is approached, 
A — > 0, signaling an instability towards a phase with 
non-zero Neel order parameter. The variation of A as 
a function of frustration is shown in Fig. 2. Let us men- 
tion that as we move close to the critical point from the 
disordered side, the density of triplets increases and is ap- 
proximately 0.3 at J2/J1 = 0.38. This leads, in principle, 
to a 30% uncertainty in the results due to the omitted 
higher order diagrams. However the accuracy in the posi- 
tion of the critical point (0.38) is much better than 30% 
because of the steep dependence of the gap on J2/J1. 
Nevertheless, within the accuracy of our calculation the 
point where A w 0.05Ji (see Fig. 2) is practically indis- 
tinguishable from the critical point. 

In Ref. [p2j it was pointed out that an additional col- 
lective low-energy mode also exists near (J2/ ' J\) c i- This 
excitation is a bound state of two triplets with total spin 
5 = 0. We proceed to investigate its properties in more 
detail. Introducing the total (Q) and the relative (q) mo- 
menta of the two triplets forming the bound state, the 
two-particle singlet is: 

l*Q> = E Q)<Q/ 2+q *U/2- q |0)' ( ? ) 

The bound state wave- function ^(q, Q) satisfies the in- 
tegral equation: 

[E S {Q) - ^q/2+q - wq/a-q] *(q, Q) = 

which can be easily derived by noticing that it is equiva- 
lent to the two-particle Schrodinger equation: H\$q) = 
E s (Q)\^q). Here ^(Q) is the energy of the collective 
mode. The function T(q) is the two-particle scattering 
amplitude from (|(j): 

rrvi \ , . (Jl + J2) , . -h , , 

1 (k) = — cos k x H cos k y + — cos k x cos k y . 

(9) 

This interaction leads to attraction between two triplets 
in the singlet channel. In Eq.(§h the (second-order) con- 
tribution of H3 , Eq. (||) into binding has been neglected. 



We have checked that the perturbative inclusion of this 
term indeed leads only to a small change of the re- 
sults presented below. Since we have to take U — > 00, 
the following replacement must be made on the right 
hand side of Eq.(|): U /dp*(p, Q) -> A, where A is 
a Lagrange multiplier, determined self-consistently from 
the condition J dq\&(q, Q) = 0. The bound state ex- 
ists only if a solution of Eq.(^) can be found such that 
E S (Q) < E C (Q) = min q [cj(Q/2+q)+w(Q/2-q)], mean- 
ing that it must be below the two-particle scattering con- 
tinuum E c . In Fig. 3 we present the numerical solution 
of Eq.(^) for a fixed value of frustration (J2/J1 = 0.4) 
above the critical value. We have also plotted the one- 
particle spectrum u>(k), and the shaded region is the two- 
particle scattering region E > E c (k). The bound state 
is stable (non-zero gap) for all k throughout the disor- 
dered phase, with a minimum of the dispersion E s (k) 
at k = (0,0). As frustration decreases the singlet gap 
E (0,0) decreases and appears to vanish at the critical 
point, as shown in Fig. 2. 

We believe that the existence of a singlet bound state 
at k = (0, 0) reflects the spontaneous nature of the dimer 
order. Indeed, we have checked that in models where 
the dimerization is explicit, i.e. due to stronger ex- 
change on certain bonds, the singlet does not exist in 
the neighborhood of k = (0,0), meaning that its bind- 
ing energy is zero. The binding energy is defined as 
e(k) = E c (k) - E s (k). We have found that e(0,0) = 
both in the 2D dimer model (as defined in Ref. ||) and 
in the two-layer model (Ref. H). Unlike the above two 
models, the Hamiltonian Eq.(0) does not break any lat- 
tice symmetries, but the ground state of Fig.l certainly 
does. It was argued in Ref. p2| that the vanishing of the 
spontaneous dimer order at the critical point is intrinsi- 
cally related to the low energy singlet mode. Thus we 
believe that the singlet bound state reflects the presence 
of non-zero dimer order parameter in the J\ — J2 model, 
i.e. the spontaneous breakdown of the discrete lattice 
symmetries in the quantum disordered phase. 

Even though we have shown that the gaps for both 
the triplet and singlet modes vanish at the critical point, 
this does not necessarily mean a departure from the 0(3) 
universality class which is related to the triplet only. In 
what follows we will in fact demonstrate that the singlet 
mode is "irrelevant" at the critical point. Let us start 
with the observation that the relevancy (or otherwise) of 
a soft excitation is directly related to its spectral weight. 
We find, as expected, that the spectral weight of the 
triplets (the residue of the one-particle Green's function) 
stays finite at the transition point. However, the spectral 
weight of the collective singlet excitation is proportional 
to the binding energy (which in turn is inversely pro- 
portional to the size of the bound-state wave function) . 
Observe that the lower edge of the two-particle contin- 
uum at k = (0,0), E c (0,0) = 2uj(Q af ) = 2A, and since 
at the transition point A — > 0, the binding energy must 
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vanish as well, e(0, 0) — > 0. The variation of the binding 
energy as a function of frustration, obtained by solving 
Eq.(||) numerically, is shown in Fig. 2. Since we realize 
that for A — > the accuracy of the calculation decreases, 
let us find the asymptotic behavior of e(0, 0) in this limit 
analytically. The one-particle dispersion around its min- 
imum has the form: w 2 (k) = A 2 + c 2 |k — Qa_f| 2 , which 
is valid close to the critical point. The triplet veloc- 
ity c is known to remain finite at the transition [(^-|l0|, 
Ac ~ Ji, where A ~ 1 (in units of the inverse lattice 
spacing) is a characteristic momentum. Denoting the 
right hand side of Eq.|| by $(q, Q), the solution of (||) 
at Q = (0,0) is: *(q,0) = $(q,0)[£ s - 2w(q)]- 1 . We 
define E s = E S {Q,Q) and e = e(0, 0) from now on. The 
energy can be found from the condition: J dq^>(q, 0) = 0. 
This integral diverges logarithmically at Qaf for small 
binding: 



d 2 q 



2a; (q) - E s 



A A, /Ac 
- + — In — 



e<A<cJi. (10) 



We remind that e = 2A — E s . When estimating the 
integral, <&(q, 0) can be replaced by <&(Qaf,0). This 
quantity is finite at the critical point which follows from 
Eq. ([)]). From (|l(J) we find the binding energy: 



e = £q exp 



i-'i) 



A <c Jx, 



(11) 



where eo> £ i ~ Ji are two constants that depend weakly 
on J2 and are finite at the critical point. The result 
is similar to the formula for the exponentially small s- 
wave bound state in a 2D potential ^6|. However there 
are two differences from the usual expression p6| : 1.) 
the pre-exponential factor in Eq.(|ll|) docs not depend 
on the "mass" A because of the relativistic form of the 
dispersion near the transition point, and 2.) because of 
the hard core constraint (U — > 00) the exponent e\ can 
not be written as | J Udr \ ~ 1 , where U is the attractive 
potential (i.e. the expression in the curly brackets in 
Eq.®). 

In Fig. 2 we present a fit of the formula (|ll|) to the 
numerical solution of (||) . The major disagreement occurs 
only close to the critical point, for J2/J1 < 0.41. The 
asymptotically exact solution predicts an exponentially 
fast vanishing of e which can not be captured reliably in 
the numerical solution. 

The size R of the bound state (in units of the lattice 
spacing) is determined by the spatial extent of the wave- 
function: R 2 = ^ q \d^f (q) / dq\ 2 , where \& is assumed 
normalized. Evaluating this expression for A — > leads 
to 
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increases as the binding energy decreases. At the crit- 
ical point R diverges exponentially with the correlation 
length. This in turn implies that the spectral weight of 
the collective singlet vanishes exponentially fast. Con- 
sequently the singlet bound state does not influence the 
triplet dynamics near the critical point and hence can 
not change the 0(3) universality class. 

Unlike the theory of Read and Sachdev p9| , our ap- 
proach does not relate directly the gap in the singlet spec- 
trum and the dimer order parameter. The latter quan- 
tity is defined as D = (S 2 .S 3 ) - (S1.S2) (see Fig.l.). We 



where £ ~ A 1 is the correlation length, and ~ 1 
is a constant. Thus we find, as expected, that the size 



have presented arguments |22| that the low-energy sin- 
glet affects the dimer order by increasing the quantum 
fluctuations. This effect becomes stronger and stronger 
as the critical point is approached and one could ulti- 
mately expect that D vanishes. High order dimer series 
results support this conclusion [p2f , but can not deter- 
mine the critical behavior of D. Thus we can not verify 
the prediction of the large- N theory (which follows from 
Eq.(Q)) that D should vanish with a large exponent. 

In summary, we have found that in the 2D J\ — J2 
model the critical behavior near the transition between 
the spontaneously dimerized and the Neeel phase is char- 
acterized by two soft modes - the usual triplet magnon 
mode, and a collective singlet excitation. Even though 
the gap in the singlet spectrum vanishes at the transition, 
we argue that it does not influence the critical dynam- 
ics of triplet excitations. The reason is that the spectral 
weight of the singlet vanishes exponentially fast at the 
quantum critical point. We have to note that our picture 
is different from that of Read and Sachdev (l^Jl^] , based 
on the large- N expansion. In jl8|,[l9| the singlet gap was 
found to be much smaller than the triplet one, which in 
turn leads to the second large length scale. In our pic- 
ture the singlet gap is approximately equal to two triplet 
gaps, and the large length scale comes from the size of 
the singlet bound sate. Nevertheless the final conclusion 
concerning the triplet critical dynamics is the same: the 
0(3) universality class describes the transition between 
the Neel ordered and quantum disordered phase. 

Finally let us mention that in cases where the ordered 
phase is characterized by incommensurate correlations 
(e.g. the triangular lattice Heisenberg model) the criti- 
cal behavior could be quite different. On semi-classical 
level the order parameter has 5*0(3) symmetry and con- 
sequently an 0(4) universality class is possible ^?J. If 
indeed the 0(4) class is realized in this case, the singlet 
mode should become a truly Goldstone mode at the tran- 
sition and therefore must have very different properties 
from the ones found in the present work. 
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FIG. 1. The Ji — J2 model on a square lattice. The cir- 
cles represent spins paired in singlets in the columnar ladder 
dimerization pattern. 

FIG. 2. Gaps in the one and two-particle excitation spec- 
tra. Solid squares, connected by a solid line, represent the 
triplet gap A = uj(Qaf), and open circles connected with a 
dashed line represent the gap of the two-particle singlet bound 
state at k = (0,0). Open squares are the singlet binding en- 
ergy, obtained by solving Eq.(|^) numerically, while the solid 
line is the fit, based on the asymptotic formula, Eq.(pd|), with 
e = 2.85Ji,ei = 1.20 Jj. 

FIG. 3. Triplet magnon excitation spectrum (solid line), 
and the singlet bound state excitation branch (long-dashed 
line) in the part of the Brillouin zone near (0,0) and Qaf- 
The shaded area represents the two-magnon scattering con- 
tinuum. All spectra are calculated diagrammatically for 
J2/J1 = 0.40. 
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